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Abstract 

A simple parton cascade model for the ultra-relativistic nucleus-nucleus collisions is proposed 
to investigate the Pauli blocking at the transition from hadronic to partonic state, i.e. QGP. A 
boost invariant study of the Pauli blocking is implemented in the Monte Carlo simulation for the 
first time. It turns out that the higher reaction energy the stronger the partonic Pauli effect is. 
It amounts to about ten percent at SPS and RHIC energies but nearly reaches twenty percent at 
LHC energy. Thus the effect is relevant to the formation and evolution of the quark-gluon plasma. 
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FIG. 1: Rapidity (left panel) and transverse momentum (right panel) distributions of vr^" in Pb+Pb 
collisions at v/s„,„=5500 GeV. 



Pauli blocking plays an important role in nucleus-nucleus collisions at low and intermedi- 



ate energies because the available phase space volume is quite small 

mm 

. However, 

in the ultra-relativistic nucleus-nucleus collisions since the available momentum phase space 
is huge and more than 95% produced particles are bosons one might despise the role of Pauli 
effect. Even if the available momentum space in ultra-relativistic nucleus- nucleus collision 
is really huge, the coordinate space, on the other hand, is strongly contracted at the earli- 
est stage of the collision. Although pions, kaons etc. do not suffer with Pauli effect, their 
precursors, u, d, and s quarks and their antiquarks, do really have all the influences of Pauli 
blocking if a Quark-Gluon Plasma (QGP) is formed and survives long enough in a limited 
region of phase space. Furthermore, since Pauli blocking is stronger for u and d quarks than 
s quark, it might benefit the production of strangeness relatively and the strong constrain 
of strange suppression could be somewhat overturned. Thus it is worth to study the effect 
of partonic Pauli blocking on the hadronic final state in ultra-relativistic nucleus-nucleus 
collisions. 

Even though our parton cascade model (see below) is crude in some aspects, the basic 
assumption is that at each nucleon-nucleon collision the produced particles (mainly pions) 
are deconfined at least for the time of the reaction, a few fm/c for instance, which is however 
very short. If a QGP is formed at the maximum compression stage we have then to test 
for the influence of Pauli principle since all the particles are fermions. In our model the 
elementary probability of particle production are taken from data 0|, thus we expect the 
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FIG. 2: Rapidity (left panel) and transverse momentum (right panel) distributions of in Pb+Pb 
collisions at v/s„,„=5500 GeV. 



model to be reasonable in describing the dynamics of particle production. However, the mean 
field effects and the dynamics of confinement are neglected and we assume that the quarks 
recombine simply when the density in phase space is small. The later, however, is a rather 
reasonable assumption since if the density is very high, the quarks are clearly deconfined. As 
we will show below, our calculations underestimate the available experimental data nearly 
a factor of 1.5, thus we expect that we are underestimating the Pauli blocking since it is 
more severe for larger densities (i.e. higher number of produced partons). More refined 
models which include the dynamics of the QGP plasma must really take the Pauli blocking 
into account. Simpler models which deal with hadrons only, should estimate the maximum 
density in phase space reached. If such density is too large then it is necessary to deal 
with the structure of the hadrons. From an experimental point of view, since the quarks 
recombine mainly in pions which are bosons, most of the Pauli blocking effect might be 
hidden. However, some important informations on the phase space density might be obtained 
from baryon-baryon coincidence measurement, especially neutrals to avoid the distortions 
of the long range Coulomb force. 

To our knowledge the quantum statistic effect in ultra-relativistic heavy ion collisions was 



really implemented in the study of bosonic enhancement of pion transverse momentum 
and was also considered formally in thepartonic transport equation S | . It should be empha- 
sized that a great effort was made in jq] after the pioneering study 9| in the field of parton 
cascade model. However, the partonic Pauli blocking effect was not really incorporated in 
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FIG. 3: Rapidity (left panel) and transverse momentum (right panel) distributions of H in Pb+Pb 
collisions at v/s„,„=5500 GeV. 
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the Monte Carlo simulation of currently active parton cascade models 
In this letter a simple parton cascade model for ultra-relativistic nucleus-nucleus collision is 
proposed investigating the effect of Pauli blocking at partonic (initialization) stage on the 
hadronic final state. The partonic Pauli blocking is implemented boost invariantly for the 
first time in Monte Carlo simulation. The amount of this Pauli effect is below ten percent 
at SPS and RHIC energies but reaching nearly twenty percent at LHC energy, due to the 
increased number of produced partons with increasing beam energy. 



Tabic 1. Particle yield in Au + Au collisions at ^sti,! — 200 and 130 GcV 
and Pb+ Pb collisions at ^3;7;r=5500 and 17.3 GcV^ 
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no feed-down correction. 

The sum of the yield of strange particles. 



The parton cascade model for ultra-relativistic heavy ion collision is generally composed 
of the parton initial state, the parton evolution, the hadronization, and the hadron evolution. 
However, the models in [l^ll^ lack the last two parts at present, for instance. The partonic 
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Pauli blocking, of course, plays a role in the former two parts only. As a first step the 
proposed simple parton cascade model is used for investigating mainly the Pauli blocking 
effect in the generation of parton initial state from the hadronic distributions at the highest 
density stage. 

, the parton initial 



There are two ways of creating the parton initial state. In 
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state was composed of partons from the mini-jets production in nucleus-nucleus collision 



S5 proauction m nucieu 
was specified in [l^ 



and the HIJING multiple mini-jet generator [15|] was specified in 
initial state in 



ll[. The parton 



12l llS^ was created via probability distributions first for the spatial and 
momentum coordinates of nucleons in colliding nuclei and then for the flavor and spatial 
and momentum coordinates of partons in nucleons. Our parton cascade model follows the 
former way, however, the JPCIAE multiple mini-jet generator |l^ is used instead of HIJING. 

The JPCIAE multiple mini-jet generator for ultra-relativistic nucleus-nucleus collision is 
based on PYTHIA [6] which is a well known event generator for hadron-hadron collision. 
In the JPCIAE model the radial position of a nucleon in colliding nucleus A (indicating the 
atomic number of this nucleus as well) is sampled randomly from Woods-Saxon distribution. 
Each nucleon is given a beam momentum in z direction and zero initial momenta in x and 
y directions (fermionic motion is neglected). The Lorentz contraction is taken into account 
after the initialization of nucleons. A least approaching distance for each colliding nucleon 
pair along their straight line trajectory is calculated together with its collision time under 



the requirement that the least approaching distance must be less than or equal to yCtot/T^- 
Here atot refers to the total cross section assumed to be equal to 50 mb for Au + Au collisions 
at a/s„„=130 and 200 GeV and 40 and 100 mb for Ph + Ph collisions at i/s„„=17.3 and 
5500 GeV 17[, respectively. The nucleon- nucleon collision with the least collision time is 
then selected from the initial collision (time) list performing the first collision. This nucleon- 
nucleon collision is modeled by PYTHIA with string fragmentation switched off, thus the 
produced particles are quark pairs, diquark pairs and gluons. To investigate the partonic 
Pauli blocking efficiently the diquark (anti-diquark) is split into quarks (antiquarks) and 
the gluons into quark pair randomly. The produced partons propagate along straight line 
trajectory similarly to the nucleons. However, the partonic interaction is neglected for the 
moment, thus partons do not collide and possibly create new partons, which is a possible 
reason for underestimating the data. After a nucleon-nucleon collision both the particle list 
and the collision list are updated and the new collision list is still composed of nucleon- 



nucleoli collisions only. The next collision is selected from the new collision list and the 
processes above are repeated until the collision list is empty. 

For each parton i among Nnew partons produced in the current nucleon -nucleon collision 
a judgment must be performed for its Pauli blocking and its unblocking probability, Punb{i), 
is calculated (see below). If the product of unblocking probabilities of Nnew produced partons 

Punb = n Pnnb{i) (l) 
i=l 

satisfies 

Punb < ^ (2) 

where ^ is a random number, the current nucleon-nucleon collision is blocked and thrown 
away. Another nucleon-nucleon collision is selected from the collision list and the processes 
above are repeated. 



In the heavy-ion co 



lision at low and intermediate energies the Pauli effect was treated 



boost uninvariantly 

mm 

. We introduce a method of boost invariant Pauli blocking 
in Monte Carlo simulation for the first time. To calculate Punb{i) we first select the partons 
with same flavor as parton i from both the Nnew partons and the particle list formed before 
current nucleon-nucleon collision. If the total number of partons picked above is denoted by 
Npick we boost all of Npiek partons to the rest frame of parton i to be consistent with the 
method of least approaching distance (this boost is not necessary indeed, we have checked 
that the discrepancy between results from with and without this boost is nearly 1%.). Two 
three dimensional cubes: one in coordinate space with size and the other in momentum 
space with size Ajj, are defined at i (as origin). The product of these two cubes, which is an 
invariant scalar I18|, is assumed to be /i^, the consequent values of and Ap are 1 fm and 
1.24 GeV/c, respectively. If the six dimensional cell spanned between partons i and j {j e 
Npick), X yij X Zij X pxj X pyj x pzj\ (invariant scalar), is within (A^ x Ap)^, i. e. the 
conditions 



< Ar/2, \yi,\ < Ar/2, \zi,\ < A,/2, (3) 



\pxj\ < Ap/2, \py,\ < Ap/2, \pzj\ < Ap/2, (4) 

are satisfied simultaneously, the occupation number of parton with flavor as parton i, 
Noccuij), is added by one. The Xij and pXj in the above equations, for instance, are the 



X component of the coordinate distance between i and j and of the momentum of j in rest 
frame oi i {i e Nnew), respectively. Thus the occupation probabihty of parton with flavor as 
parton i reads 



where g=6 is the spin and color degeneracies of the quark and/or antiquark with a given 
flavor (identified as u, d, and s). The unblocking probability of parton i is then 



The above partonic system is then hadronized by JETSET [6| after the nucleon-nucleon 
collision ceased (freeze-out) and the hadronic final state is the consequence. One investigates 
the partonic Pauli blocking in the generation of the initial parton state and consequently on 
the hadronic final state by comparing the calculations with partonic Pauli blocking to those 
without. 

In Tab. 1 is given the multiplicity of particles in Au + Au collisions at ^s„„=130 and 
200 GeV (5% most central) and Pb + Pb collisions at 17.3 GeV (5% most central) and 5500 
GeV (10% most central) from calculations with and without Pauli blocking. The decreasing 
percentage of the particle multiplicity in the calculations with Pauli blocking relative to the 
ones without Pauli blocking is given in the table as well. In the last row the strangeness 
multiplicity, the sum of all the strange particles in the table, is shown. One sees in Tab. 
1 that the partonic Pauli blocking effect is generally below ten percent at SPS and RHIC 
energies but reaching nearly twenty percent at LHC energy. The higher reaction energy 
the stronger Pauli effect because of the increasing number of produced quarks and of the 
competition between the partonic Pauli blocking, which decreases pion production, and the 
fragmentation, which benefits low pt pions. The partonic Pauli effect is also of benefit to 
the production of strangeness as seen in Pb + Pb collision at LHC energy. 

The rapidity (left panel) and transverse momentum (right panel) distributions with and 
without Pauli blocking for Pb + Pb collisions at y/Snn= 5500 GeV are given, respectively, in 
Fig. in El andElfor n'^, and S~, for instance. The discrepancy between the two cases is 
visible, in rapidity distribution especially. 

In summary, a simple parton cascade model for ultra-relativistic nucleus -nucleus colli- 
sions is proposed studying the partonic Pauli blocking effect in the parton initial state and 
consequent hadronic final state. A boost invariant study of the Pauli effect is implemented 



Poccu{i) = Noccu{i)/g 



(5) 



Punb{i) = 1 -Voccuii)- 



(6) 
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in the Monte Carlo simulation for the first time. It turned out that the amount of this Pauli 
blocking effect is below ten percent at SPS and RHIC energies but reaching nearly twenty 
percent at LHC energy, i.e. the higher beam energy the stronger the partonic Pauli effect 
is. Thus the effect is relevant to the formation and evolution of the quark-gluon plasma. 
However, the results from our model are a kind of under estimation. Since the charge mul- 
tiplicity in Au + Au collisions at ^s„„=200 and 130 GeV in the calculations with Pauli 
blocking is, respectively, 3063 and 2562 (including feed-down correction from k^, A, and A) 
and the multiplicity of 7r+ and vr" in Pb + Pb collision at i/s„„=17.3 GeV is, respectively, 
416.5 and 429.3 in Tab. 1. They are all less than the corresponding experimental data 
of 4630±370 Q, 3860±300 Q, and 619±17±31 and 639±17±31 0, respectively. On 
the other hand, the results of our parton cascade model are not too bad to compare with 
experimental data, which means that the model is reasonable in its range. 
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performed) are acknowledged. 
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